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We present Raman measurements on the iron-pnictide superconductors CeFeAsOi-xFcc and 
NdFeAsOi-ccFx. Modeling the Fe-As plane in terms of harmonic and a cubic anharmonic Fe-As 
interaction we calculate the temperature dependence of the energy and lifetime of the Raman active 
Fe Big mode and fit to the observed energy shift. The shifts and lifetimes are in good agreement 
with those measured also in other Raman studies which demonstrate that the phonon spectrum, at 
least at small wave numbers, is well represented by phonon-phonon interactions without any signifi- 
cant electronic contribution. Even at zero temperature there is a non negligent effect of interactions 
on the phonon energy, which for the Fe Big mode corresponds to 6cm~^, or 3% of the total energy 
of the mode. We also estimate the anharmonic expansion from Fe (56 54) isotope substitution 
to Aa ^ 5.1 • 10~^A and Adpe-As ~ 2.5 • 10~^A and the shift of harmonic zero point fluctuations 
of bond lengths < Ax^ >< 3 • 10"^ A^ giving a total relative average decrease of electronic hopping 
integrals of \6t\/t < 2.0 • 10"^. For a non-phonon mediated weak coupling superconductor this gives 
an isotope exponent a ^ 10~^. The results poses a serious challenge for any theory of supercon- 
ductivity in the pnictides that does not include electron-phonon interactions to produce a sizable 
Fe- isotope effect. 



I. INTRODUCTION 

The recently discovered iron-pnictide high tempera- 
ture superconductor^^ have several features in common 
with the much studied cuprate superconductors. In 
addition to their overall quasi-two dimensional nature 
such features include commensurate antiferromagnetism 
in close proximity to or even coexisting with supercon- 
ductivity. Intriguingly, there is also a common tendency 
for breaking of the fourfold rotational symmetry of the 
crystal^ which in the case of the cuprates is most likely 
an intrinsic electronic property related to the strong 
correlations.'^ For the pnictides (as well as in some cases 
for the cuprates) there is at low doping a crystal sym- 
metry breaking from tetragonal to orthorombic which is 
closely connected to an accompanying spin density wave 
(SDW) order. ^ ^ This interplay between atomic and elec- 
tronic degrees of freedom naturally leads to the question: 
What is the role of phonons for the electronic proper- 
ties of these systems? Focusing on the pnictides, re- 
sults from density functional theory (DFT)^ ^ are cred- 
ible as they give an electronic spectral distribution con- 
sistent with photoemission experiments^ but have also 
found that the electron-phonon coupling is much too 
weak to explain the high TcJ Challenging these find- 
ings was the recent report of an isotope shift of Tc (and 
T^SDw) when substituting ^^Fe with ^^Fe with an expo- 
nent of a = -4t^ = 0.4 close to the ECS value of 0.5 

a mm 

for a pure iron mode.!^ The mechanism for superconduc- 
tivity in the pnictides is a major unresolved issue and it is 
clearly essential to understand better the lattice-electron 
interplay. 

In this paper we measure and model one particular Ra- 



man active phonon which lives primarily on Fe atoms and 
thus common to all iron-pnictide superconductors. The 
purpose is to study the anharmonic structure of the Fe- 
As plane with two main objectives: First, to investigate 
whether any non-phonon contributions are necessary to 
describe the temperature dependence of the phonon en- 
ergy and lifetime. Second, to estimate the magnitude 
of the lattice expansion that follows from substitution 
with a lighter mass and the corresponding changes in 
electronic hopping integrals, thus investigating the plau- 
sibility o f a n on-phonon related isotope effect in the 
pnictides.ESEH 

In brief, we model the harmonic spectrum of the 
As-Fe plane and use standard Greens function meth- 
ods to calculate the self energy of the Fe Bi^ Raman 
phonon to second order in the cubic anharmonic Fe- 
As interatomic coupling. The harmonic spectrum is 
based on a minimal parameter fit to spectra derived 
within the local density approximation (LDA) of DFT 
but the calculation of the temperature dependent contri- 
bution from phonon-phonon interactions is beyond the 
capabilities of that method. We present Raman re- 
sults on (Nd,Ce)FeAsOi 

—x^x ^iid estimate the anhar- 
monic coupling strength by fitting our model calculations 
to the measured temperature dependent energy shift. 
Also the width agrees well with Raman measurements 
on CaFe2As2^^ thus effectively ruling out any significant 
electronic contribution to the broadening. Based on the 
magnitude of the anharmonic coupling we estimate the 
isotope (^^Fe^^^Fe) shift of the lattice parameters to 
< 2 • 10~^ and calculate a similar relative decrease of 
interatomic hopping integrals. Harmonic zero point fluc- 
tuations give an even smaller isotope shift of hopping 
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ble subtractive mode. In all scans the 514.5 nm line from 
a Ar+ laser used with a power of less than ImW was 
focused onto the samples with a spot size less than 2 jiia. 
The samples were installed in a LHe cooled cryostat. We 
observe three Raman active modes with energies around 
170, 210 and 220 cm-\ see Fig. fljand Fig. ]^ in agree- 
ment with other Raman studies, that are identified as 
(Ce/Nd)-Ai^, As-Ai^, and Fe-Bi^ modes respectively. 
Also in agreement with earlier studies we found no effect 
on the phonon energies from crossing into the supercon- 
ducting phase (Tc =35K/45K for Ce/Nd). As discussed 
in the introduction we will be interested in analysing the 
temperature dependence of the energy as extracted in 
Fig. [3] and for modeling purposes we consider only the 
Fe-Bi^ {x^ — y'^) mode which lives primarily in the Fe-As 
plane. 



FIG. 1: (Color online) Temperature dependent Raman spec- 
tra for NdFeAsOi-ccFcc (x = .12) showing three Raman active 
phonon modes labeled by their symmetry and main atomic 
displacement. 



integrals. In weak coupling theory correspondingly small 
changes in the electronic density of states (DOS) is ex- 
pected to give an isotope exponent a ~ 10~^. 



II. RAMAN SPECTRA 

Raman spectra between lOOcm"^ and 400cm~^ were 
collected at temperatures ranging from 20K to 300K for 
polycrystalline samples of CeFeAsOi_a^Fa^ [x = 0.16) 
and NdFeAsOi_a,Fa, {x = 0.12).^ For sample prepara- 
tion and characterization see Chen at al.^^ All spectra 
were recorded using a Dilor-XY800 spectrometer in dou- 
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FIG. 2: (Color online) Raman spectra for CeFeAsOi-xF^ 
{x .16). 
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FIG. 3: (Color online) Temperature dependence of the 
phonon peak position for NdFeAsO0.88F0.12 (red) as extracted 
from Fig. [lland correspondingly for CeFeAsO0.84F0.i6 (blue) 
from Fig. [2] 



III. MODELING 

The temperature dependence of the lifetime and en- 
ergy of a phonon is ordinarily due to phonon-phonon in- 
teractions that arise from anharmonic interatomic poten- 
tials. In general the cubic anharmonicity is the dominant 
term^^ and we consider only this. The Raman {q = 0) in- 
tensity for Stokes scattering for the mode j at frequency 
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uj is given by Is{j,^) oc -{l^n{uj))ImVret{q = 0, 
where Vret is the retarded phonon Greens function which 
to hnear order in the self energy H = A — iT (minus sign 
by convention) gives 

r t- (i + nH)r(o,j» 

1S\1-,0J) CX =5^ ^ =5^ (1) 

[u;o(0,i) + A(0,i,c.)-u;]2 + r2(0,i,c.) 

where n{ijj) = (e^^/^^^ — 1)"-^ is the Bose occupation 
factor. The measured width (FWHM) is thus ideahy 



given by V and the shift by A. We calculate the 
self energy to second order in the interaction Ha = 

where A^j = a^j + ^-(f j phonon operator.!^ V 

are the matrix elements 



ri,r2,r3,Ri,R2,R3 



n ( 

.1=1 ^ 



2Mr,uj{qiJi) 



2 p'^Qi '^i p'-y 



)2e 



(?'l,?'2,r3,-Rl,i?2,-R3) 



(2) 



where the derivatives act on ^, the interatomic poten- 
tial, with Ti the atomic positions within unit cell Ri and 
^qiJii'^i) are the displacement vectors of the respective 
mode. The expression thus amounts to the third derivate 
with respect to phonon displacements of the atomic po- 
sitions. The iron atoms are only weakly coupled to the 
{Ce^Nd)0 layer and since the phonon dispersions in 
the c-direction are wealP we can model the the Fe-Bi^ 
mode by considering only the Fe-As plane. (This also 
makes the procedure quasi-universal as it only depends 
on the properties of this plane which is common to all 
iron-pnictide superconductors.) For the harmonic prob- 
lem we use a minimal spring model (Fig. [4| with near- 
est neighbor Fe-As potential |(5r^, nearest neighbor Fe- 

Fe potential ^5r''^^ and in-plane As- As coupling k^8r"'^ 
with (5r, 5r' ^ and 5r" the deviations from the respec- 
tive equilibrium distances. For the Fe-As coupling we 
also add a cubic term — |(5r^, where the magnitude of g 
is to be estimated from the experimental fit. The Ra- 
man active Bi^ is here a pure Fe mode as depicted in 

Fig. [4] with energy us^g = \[^^^ = 56^^) and 
to get ujB^g = 220cm~^ we take k = 8.7eV/A^, us- 
ing = 35^ W Values k' = 0.3k and k" = 0.2k gives a 
lower edge of the accoustic branches in good agrement 
with LDA calculations.^ Using this phonon spectrum 
we calculate the interaction matrix elements through 
Eq. [2] and the cubic potential. Sampling the coupling 
V{0,FeBig]qJi] -qj2) = Vo{qJij2) over the Brillouin 
zone shows that except at high symmetry points all 12 
modes are coupled quite isotropically and does not vary 
strongly with q if the frequency dependence is included, 
giving 

^o(^ Jl, J2) ^ ^g{-^f^' ^ Vo , (3) 
^mcjB^g 



with the numerical value k, = 0.10. 




FIG. 4: (Color online) Schematic of the Fe-As plane with 
in-plane lattice parameter a. The harmonic interatomic cou- 
plings are /c, and k'' and cubic coupling g as described in 
the text. Arrows indicate the Fe Bi^ mode. 



A. Phonon self energy 



Now we are ready to address the frequency shift by 
calculating the self-energy to second order in Vq. The self 
energy is given by the bubble diagram^^ with imaginary 
part 
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FIG. 5: (Color online) Phonon density of states of the isolated 
Fe-As plane calculated within the harmonic model discussed 
in the text and its approximation in terms of three Lorentzians 
(solid curve) that contain 4,3,5 modes with increasing energy 
respectively. 

(The phonon-difference process, the second term, van- 
ishes at zero temperature but can give a significant con- 
tribution at finite temperature.) The simplest standard 
way to evaluate this expression is to assume that the 
scattering is diagonal in the modes Vq ~ ^jus^ 
mens model, which gives a characteristic temperature 



The expression is naturally understood as the scatter- 
ing of the mode at into two modes within the same or 
different DOS peaks (or the corresponding difference pro- 
cess). Note that ro(cc;Big) will be dominated by scatter- 
ing into the low-energy peak (i = j = 1) which will give 



dependence r(cj) ~ (1 + 2n(cc;/2)). As discussed previ- 
ously we find instead that the scattering is approximately 
isotropic between modes and we need to do a more care- 
ful analysis. In principle we could calculate this expres- 
sion numerically at each temperature using our numerical 
phonon spectrum but instead we will use an approximate 
Lorentzian fit for the phonon DOS with the advantage of 
giving an analytic expression for both F and A including 
the full temperature dependence. The calculated phonon 
DOS, plotted in Fig. |5) consists of three regions of high 
density centered around wi ~ lOOcm"^, W2 ^ 200cm~^, 
^ 300cm~^. The lower part contains mi = 4 modes, 
the intermediate part contains m2 = 3 modes and the 
upper part ma = 5 modes giving 

TT (W — Wi)'^ -\- "y"^ ^-^ 

i=l i 

which is thus normalized by p{ijS) — 12N and where 
we estimate 7 = 25cm~^. To proceed we evaluate the oc- 
cupation factor n{u;) at the peak position of the respec- 
tive DOS and assume a q-independent spectral weight 
giving F as an integral over all pairs of peaks of the spec- 
tral weight with the result 



(6) 



a temperature dependence close to the Klemens model 
F(a;) ~ (1 + 2n(a;/2)).CSl Using the Kramers-Kronig rela- 
tion A(lj) = — - ^^dcj' gives 

^ ^ TT J— 00 UJ' —LJ ' " 
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+ (n(cji) - n{uj))[ 
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FIG. 6: (Color online) Fit of the model frequency shift (Ao) 
to data for NdFeAsO.ssF, 12 (boxes) and CeFeAsO .qaF.iq 
(diamonds). Inset are the corresponding model lifetimes (Fo) 
together with linewidth data (triangles) on CaFe2As2 from 
Choi et al.l^. 



IV. DISCUSSION 

Figure [6] shows the fit from equation [7| to the exper- 
imental values using a; = cjq + Ao(a;Big), with fitting 
parameters uoq and g (Aq ^ g'^)- For the Nd sample we 
find uoq = 221. 7cm-^ and g = 103eV/A^ and for the Ce 
sample coq = 221.0cm~^ and g = 116eV/A^, where uoq is 
close to the assumed harmonic value 00 Big = 220cm~^^ 
It is interesting to note that we find Ao(T = 0) ^ 6cm~^, 
i.e. even at zero temperature phonon-phonon interac- 
tions give a finite (here 3%) contribution to the phonon 
energy, an effect which LDA calculations of phonon en- 
ergies based on linear response neglects. 

Due to the polycrystalline nature of our samples we 
are not able to find reliable phonon widths and compare 
instead our theoretical results to Raman studies by Choi 
et al^ on CaFe2As2. In that material the Bi^ mode has 
lower energy but the temperature dependence is similar 
and the width correspond well with our calculations with- 
out any additional fit. The temperature dependence is in 
fact just the Klemens model but the strength of the cal- 
culation is that only one parameter g gives both the shift 
and the width. In Rei'^^J the linewidth variations were 
tentatively assigned to changes in the electronic scatter- 
ing below TsDW' These speculations we can quite defi- 
nitely rule out; the anharmonic contribution accounts to 
good accuracy for both the shift and linewidth variations 
of this mode.— 



A. Lattice expansion 

The cubic coupling will give rise to an expansion of 
the lattice that we estimate by considering an isolated 



Fe-As bond for which < 0\6r\0 >- 4(^^^1/2^3/2 
m* = Mm/{m + M). We include also a nearest neigh- 
bor Fe-Fe anharmonic coupling g' with the same relative 
strength {g' = g{k' /k)^^'^). With this the isotope sub- 
stitution ^^Fe^ ^^Fe gives an expansion of the in-plane 
lattice parameter Aa = 5.1 • 10""^ A and Fe-As distance 
Adpe-As = 2.5 • 10~^A, which compared to a ~ 4A and 
dpe-As ^ 2. 5A.'^ gives a relative expansion < 2 • 10~^. 

What is the possible effect of a small lattice expan- 
sion on Tc and Tsdw^ There is a direct effect on the 
electronic hopping integrals and theories of an iso- 
tope effect based on this has been suggested for the 



cuprates and Ceo as well the pnictides.^^^^I^ Assuming 
t — ^06"^*^^/^°"-^^ with q ~ ll^and where f = ror + (5 with 



small displacement S we find = 



2 



2 r2 



where 6r = S • r and 5± = 6 x r. The linear term has 
only an anharmonic contribution whereas the quadratic 
terms will get contributions from zero point fluctuations 
in the harmonic approximation. We estimate 5'^ from 
the ground state energy per atom eo = 38.9meV and 
eo = 39.2meV for ^^Fe and ^^Fe respectively and energy 
eo/2 per Fe-As bond. The fluctuation along a bond is 
given by |(^r^ ^ eo/4 {k = 8.7eV/A^), giving the differ- 
ence (5^ ^ Acq /{2k) ^ 2.3- 10~^A^ and for the transverse 
fluctuations 5l ^ 28l ^ 4.6 • IQ-^A^^s rj.^^ contribution 
to the shift of the hopping integrals is partially canceled 
by the different signs and in total smaller by an order of 
magnitude compared to the anharmonic contribution. 



B. Isotope effect 

Consider a weak coupling (non phonon mediated) SC 
or SDW transition with Tc - Qe'^^^^^^^ , where n is 
the relevant energy cut-off, which may be the magnon 
energy for spin fluctuation mediated pairing, V is the 
effective interaction and N{0) the relevant electronic 
DOS at the Fermi energy. Focusing on the contribution 
from the DOS, A^(0) ~ 1/t gives the isotope exponent 

^ dlnTe _ 1 dlnNjO) _ O.S-IQ-^ _ -|r,-2 /Ao 

" ~ dlnm ^ N{0)V dlnm ^ N{0)V ^ ' 

suming, N{0)V ^ 0.5. It cannot be much smaller to get 
a high Tc.) The interaction strength and Q may also de- 
pend on the hopping integrals and give isotope shifts of 
similar magnitude. Although this is a simple analysis we 
expect the order of magnitude estimate to be relevant to 
any purely electronic microscopic model containing inter 
and intra orbital hopping integrals and interactions that 
only depend indirectly on the lattice parameters. 

Alternatively, we may relate the change of lattice pa- 
rameter to a corresponding pressure of dP = ^^/P ^ 
0.4kbar through the compressibility (3 = — ^^p^ ^ 1.0 • 
10~^/kbai'^ (We have no estimate of the c-axis change 
related to the in-plane expansion but only assume that 
this is of similar relative magnitude.) Pressure depen- 
dence of Tc of around 0.2i^r/kbar has been reported in 
several materials at different dopings although close to 
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optimal doping it appears that the effect is generahy sig- 
nificantly smaller. Nevertheless, from these considera- 
tions we find an upper estimate of the isotope exponent 
a ^ 0.06. 

Clearly something more sophisticated is needed to pro- 
duce a = 0.4 as found experimentally in RefEl Encour- 
aging perhaps, the sign from the naive weak coupling 
analysis based on a change in the electronic DOS does 
agree with experiments and would naturally imply a sim- 
ilar exponent for both SC and SDW.-^ 

V. SUMMARY 

In summary, we find that the temperature dependent 
shift and width of a Raman active phonon is well repre- 
sented by the anharmonic contribution, consistent with 
weak electron-phonon interactions. At the same time, we 
estimate the change of electronic hopping integrals due to 
Fe isotope substitution and find that these are too small 
to generate an isotope effect on Tc or Tsdw of the mag- 



nitude reported without an explicit phonon contri- 
bution. These apparently contradictory results present a 
significant challenge for any theory of superconductivity 
in the iron-pnictides. 
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perimental error of ~ 1 • lO^^A.I^ 
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